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ABSTRACT 


We have observed the z = 4.3 protocluster SPT2349—56 with the Australia Telescope Compact Array 
(ATCA) with the aim of detecting radio-loud active galactic nuclei (AGN) amongst the — 30 submil- 
limeter (submm) galaxies (SMGs) identified in the structure. We detect the central complex of submm 
sources at 2.2 GHz with a luminosity of Lə > = (4.424+0.56) x 10% WHz-!. The Australian Square 
Kilometre Array Pathfinder (ASKAP) also detects the source at 888 MHz, constraining the radio spec- 
tral index to a= —1.6-£0.3, consistent with ATCA non-detections at 5.5 and 9 GHz, and implying 
Li 4, vest = (2.40.3) x 1026 W Hz ^!. This radio luminosity is about 100 times higher than expected 
from star formation, assuming the usual far-infrared (FIR)-radio correlation, which is a clear indica- 
tion of an AGN driven by a forming brightest cluster galaxy (BCG). None of the SMGs in SPT2349—56 
show signs of AGN in any other diagnostics available to us (notably '*CO out to J — 16, OHi3um, 
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[C t] /IR, and optical spectra), highlighting the radio continuum as a powerful probe of obscured AGN 
in high-z protoclusters. No other significant radio detections are found amongst the cluster members, 
with stacking on either all members or just the ten most luminous members yielding non-detections 
consistent with the FIR-radio correlation for star-forming galaxies. We compare these results to field 
samples of radio sources and SMGs, along with the 22 SPT-SMG gravitational lenses also observed 
in the ATCA program, as well as powerful radio galaxies at high redshifts. Our results allow us to 
better understand the effects of this gas-rich, overdense environment on early supermassive black hole 
(SMBH) growth and cluster feedback. We estimate that (3.3-- 0.7) x 10% W of power are injected into 
the growing intra-cluster medium (ICM) by the radio-loud AGN, whose energy over 100 Myr is com- 
parable to the binding energy of the gas mass of the central halo. The AGN power is also comparable 
to the instantaneous energy injection from supernova feedback from the 23 catalogued SMGs in the 
core region of 120 kpc projected radius. The SPT2349—56 radio-loud AGN may be providing strong 


feedback on a nascent ICM. 


Keywords: Submillimeter astronomy (1647) — Galaxy evolution (594) 


1. INTRODUCTION 


Submillimeter (submm) galaxies (SMGs) are impor- 
tant sites of stellar mass build-up at cosmic noon and 
earlier (e.g., Chapman et al. 2003, 2005; Smail et al. 
2004), with star-formation rates (SFRs) as high as hun- 
dreds to thousands of solar masses per year. Several 
studies have also suggested that SMGs may be good 
tracers of dark matter halos at early cosmic time (e.g., 
Blain et al. 2004; Chen et al. 2016; Dudzeviciute et al. 
2020). Simulations conducted by Miller et al. (2015) 
found that while many dark matter halos at z — 2-4 do 
not contain any SMGs, large and rare associations of five 
or more SMGs do trace massive overdensities of dark 
matter that have the potential of evolving into present- 
day massive clusters. Supporting this, in the recent past, 
several high-redshift protoclusters have been identified 
entirely through their submm emission (e.g., Chapman 
et al. 2009; Daddi et al. 2009; Capak et al. 2011; Casey 
et al. 2015; Miller et al. 2018; Oteo et al. 2018; Gómez- 
Guijarro et al. 2019; Wang et al. 2021). 

AGN and star-formation processes in galaxy evolution 
are clearly related (e.g. Kormendy & Ho 2013). En- 
hanced AGN activity relative to the field environment 
has been found in massive protoclusters at z = 2-3 (e.g., 
Pentericci et al. 2002; Lehmer et al. 2009; Digby-North 
et al. 2010), which is likely related to the enhancement 
of star formation in galaxy protocluster members (e.g. 
Elbaz et al. 2007; Chapman et al. 2009; Brodwin et al. 
2013; Casey et al. 2015; Gilli et al. 2019). The sup- 
pression of star formation in galaxy clusters requires 
mechanical and radiative feedback, which is naturally 
generated by AGN. Extended X-ray emission has been 
detected in clusters, showing empty regions or cavities 
in the hot gas (e.g. Fabian 2012), which can naturally 
be explained as shocked gas from the feedback. The 


Clusters Around Radio-Loud AGN (CARLA) survey of 
around 400 high-redshift radio galaxies (HzRGs) from 
z = 1-3 (Wylezalek et al. 2013) showed that in the ma- 
jority of cases, the radio AGN is located near the cen- 
ter of the galaxy overdensity as traced by their stellar 
mass (Spitzer-IRAC emission). This is strong evidence 
that radio galaxy feedback in a growing brighest cluster 
galaxy (BCG) is important for the evolution of massive 
galaxy clusters. 

Galaxy overdensities in the high redshift Universe 
have likely not yet virialized (e.g., Overzier et al. 2013). 
'They have abundant reservoirs of cold gas to supply star 
formation, while the ongoing mergers between galax- 
ies expected in the dense environments provide trig- 
gers for star formation. Mergers can also provide the 
tidal torques necessary for the gas to overcome its an- 
gular momentum and fall to accretion disk of the su- 
permassive black hole (SMBH). AGN require this nu- 
clear accretion as a power source. This is in contrast to 
low redshifts, where structures are virialized, and both 
AGN and star formation are largely suppressed in clus- 
ter galaxies (Ehlert et al. 2014; Rasmussen et al. 2012; 
van Breukelen et al. 2009; Martini et al. 2006; Kauff- 
mann et al. 2004). Studies of AGN in protoclusters has 
recently become a viable endeavor, with relatively deep 
Chandra and XMM- Newton observations at z — 1.5-3 
(e.g., Digby-North et al. 2010; Wang et al. 2013; Travas- 
cio et al. 2020). Continuing to study the rich variety 
of protoclusters and extending these studies to earlier 
times can inform how host galaxies are affected by their 
SMBHs, as well as the connection to the surrounding 
environment. In the z — 3.09 SSA22 protocluster, 5096 
of the SMGs were found to host X-ray luminous AGN 
(Umehata et al. 2019) — a clear excess over the 15% 
found for field SMGs (e.g. Wang et al. 2013). At larger 
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distances, an overdensity of ten SMGs found by the Her- 
shel Space Telescope at z=4.0 (Oteo et al. 2018) has 
been studied by Chandra in the X-ray (Vito et al. 2020) 
and in the radio (Oteo et al. 2018), revealing no sig- 
nificant excess of AGN activity in the system over field 
SMGs (22% versus 15%, respectively). 

The 2,500 deg? survey conducted by the South Pole 
Telescope (SPT — Vieira et al. 2010; Everett et al. 
2020) at 3.0mm, 2.0mm and 1.4mm has uncovered 
a small population of nine millimeter sources ranging 
from z=3-7, which are extremely luminous, yet ap- 
parently not gravitationally lensed (e.g. Spilker et al. 
2016; Reuter et al. 2020; Wang et al. 2021). A well 
characterized example of this is SPT2349—56, a proto- 
cluster system at z —4.303 (Miller et al. 2018). Ob- 
servations at 870 um using the Large APEX BOlome- 
ter CAmera (LABOCA ; Siringo et al. 2009) on the At- 
acama Pathfinder Experiment (APEX; Giisten et al. 
2006) telescope (with a 19-arcsec beam size) first re- 
vealed an extended structure with two distinct lobes 
connected by a bridge with a combined flux density 
of Ss70 um = (106 + 8) mJy (Miller et al. 2018; Wang 
et al. 2021). Follow-up observations with the Atacama 
Large Millimeter-submillimeter Array (ALMA; Woot- 
ten & Thompson 2009) measured the redshift of its 
brightest central source through !?CO lines (Strandet 
et al. 2016), and then resolved the structure into over 30 
submm-luminous sources (Miller et al. 2018; Hill et al. 
2020; Rotermund et al. 2021), with a velocity disper- 
sion suggesting a central halo mass of around 1013 Mo. 
A VLT/MUSE observation reveals the presence of a 
Lya blob (LAB), with a linear size of about 60kpc, 
close to the core of SPT2349—56 (Y. Apostolovski et 
al. in prep.). None of the other protocluster SMGs were 
detected as Lya emitters (LAEs) in the MUSE data. 
Lya halos are commonplace in most HzRG protoclusters 
(Venemans et al. 2007). Similar objects are often found 
in protoclusters identified through other means, for ex- 
ample optical galaxy overdensities (Overzier 2016), and 
indicate the presence of significant amounts of neutral 
gas in the assembling cluster. 

'This paper presents a search for radio detections of 
members of the SPT2349—56 cluster. Section 2 de- 
scribes the radio and (sub)millimeter observations. Sec- 
tion 3 presents the results derived from the source ex- 
traction and analysis. In Section 4 we discuss the de- 
tected central radio source, the energy injected into 
a growing ICM, and the implications for radio-loud 
AGN in protoclusters. We conclude in Section 5. 
Throughout our analysis, a Hubble constant of Ho = 70 
kms! Mpc-! and density parameters of Q4 =0.7 and 


Table 1. ATCA-selected sources within a radius of 1 Mpc 
in projection (140") of SPT2349—56. 


ID RA Dec Freq. Flux 
(GHz) (uJy) 

ID1 = = 8.98 < 159i! 

ID1 = = 5.47 < 12017 


ID1 23:49:42.760 —56:38:25.05 2.17 214027 
ID1! 23:49:42.55 —56:38:19.4 0.888 867189 
ID2 23:49:38.838 —56:37:09.63 2.17 547+36 
ID2 23:49:38.750 —56:37:06.09 0.888 1324+182 
ID3 23:49:43.692 —56:38:01.82 2.17 13526 


+ ASKAP measurement 
tf 3e ATCA limit 


Table 2. ALMA observing programs used for follow-up 
analysis. Details on additional ALMA Band-7 observations 
used in this paper can be found in (Hill et al. 2020). Here 
the frequency is the central frequency between the upper and 
lower sidebands, the continuum sensitivity is calculated at 
the center of the primary beam and averaged over the upper 
and lower sidebands, and the beam is the average circular 
synthesized beam FWHM. 


ID Date Freq. Octm Beam Array 
(GHz) (uly) (7) 


2015.1.01543.T 03/20/16 148.3 10 0.88 | C36-2/3 
2018.1.00058.5 10/03/18 146.8 12 0.28 C43-6 
2021.1.01313.5 07/27/22 146.3 21 0.27 C-6 
2021.1.01313.8 09/01/22 231.9 31 0.47 C-4 


Om —0.3 are assumed, resulting in a proper angular 
scale of 6.88 kpc/" at 2 = 4.3. 


2. DATA 
2.1. ATCA observations 


SPT2349—56 was observed by the Australia Tele- 
scope Compact Array (ATCA) at 2.2, 5.5, and 9.0 GHz 
between 2012, January 23 to 27, as part of a pro- 
gram (C1563) to observe 23 SPT-SMGs (described in 
Appendix C). We used the Compact Array Broad- 
band Backend (CABB) configured in the 1M-0.5k mode, 
which leads to a bandwidth of 2GHz per correlator 
window with 1 MHz per channel of spectral resolution. 
'The observations were performed in the most extended 
ATCA configuration, 6A, with six working 22 m anten- 
nas. The on source time was 34min, which was typ- 
ical for all SPT-SMGs observed (see Table 5). The 
data were edited, calibrated, and imaged using the 


Dec 


4 CHAPMAN ET AL. 


SPT2349 — 56 
ATCA 2.2 GHz 


-56:37:00 


30 


30 


39:00 


a 
23:49:45 40 35 30 
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Figure 1. Background: ATCA 2.2 GHz imaging of the SPT2349—56 region, with gold contours highlighting the 106 mJy 
extended LABOCA source at 870um. The linear ATCA feature east of SPT2349—56 is from the synthesized beam structure 
of a bright 20 mJy source to the south (see Appendix A). An ATCA radio source is identified near the LABOCA core, with an 
ASKAP source (white contours) overlapping. The bright ATCA+ASKAP source to the northwest is identified with a Milky 
Way star (ID2 in Table 1). Inset: A 20" x 20" zoom-in of ALMA 350 GHz continuum imaging (Hill et al. 2022) with overlays 
of ATCA 2.2 GHz (cyan), ASKAP 888 MHz (white). ATCA contours start at 3.70 revealing the FWHM of the source (4" x 8"). 
White contours (ASKAP) start at 3c, and the FWHM of the source is 16" x 25". ALMA sources are named from Miller et al. 
(2018) in order of their 850 jum flux density. The ATCA radio detection of the B-C-G complex of galaxies is evident. 
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Figure 2. Spectral energy distribution showing the ATCA 
and ASKAP radio detections at rest wavelengths, and the 
flux densities of the brightest two of the three central SMGs 
(B - red circles; C — black circles). Also shown is the rest- 
frame optical photometry of ALMA source C (as in Fig. 1), 
which was modelled with a 3 x 10! Mọ stellar mass fit 
(Rotermund et al. 2021). Source B is undetected at these 
wavelengths. The Arp220 SED (dashed blue line) is normal- 
ized to the submm photometry, revealing that the SPT2349 
BCG galaxy complex has significant excess in radio above the 
far-infrared (FIR)-radio correlation for star-forming galax- 
ies. The 5.5 GHz and 9 GHz limits are shown at 30. The 
radio spectral index is constrained to a= — 1.58 + 0.31 (fit- 
ted line), primarily by the ASKAP detection, and consistent 
with the upper limits. The grey shadings show the 1 and 20 
uncertainties in the fit (Appendix C); an Arp220 a= — 0.8 
spectral index is ruled out at the 30 level by the 5.5 GHz 
non-detection. 


Miriad package. Data affected by known radio fre- 
quency interference (RFI) or with bad visibility ranges 
were flagged accordingly. We estimate an absolute cal- 
ibration uncertainty of 596 at 2.2 and 5.5 GHz, and 
1096 at 9.0 GHz. We inverted the visibilities using nat- 
ural weighting, leading to beam sizes of 7.7" x 4.2", 
3.2" x 2.1", and 2.0" x 1.3" at 2.2, 5.5, and 9.0 GHz, 
respectively, with associated RMS noise values of 27, 
40, and 53 uJy beam !, respectively. Figure 1 displays 
the ATCA 2.2 GHz map surrounding SPT2349—56, re- 
vealing a well-detected (8e) source near the core of 
SPT2349—56. No sources at 5.5 or 9.0 GHz are found in 
the vicinity of SPT2349—56. ATCA sources surround- 
ing SPT2349—56 out to 1 Mpc in projection are listed 
in Table 1, and the wider-field ATCA map is shown in 
Appendix A. 


'The shortest baseline is 30m and the images should 
be sensitive to emission on angular scales up to a few 
arcminutes. In principle, these data should not be miss- 
ing any flux on the scales covering both the ATCA and 
ASKAP (see below) sources, although the ATCA data 
will be less sensitive to lower surface brightness emis- 
sion. However the short 34 min integration, with quite 
limited SNR, may still be missing some structure due to 
sparse uv coverage. Similar issues were discussed in an 
ATCA snapshot survey of distant HzRGs (De Breuck 
et al. 2000) and are elaborated in section 4. 


2.2. ASKAP observations 


The Australian Square Kilometre Array Pathfinder 
(ASKAP) comprises 36 twelve-metre dishes located in 
the Inyarrimanha Ilgari Bundra! at the CSIRO Murchi- 
son Radio-astronomy Observatory (MRO) in Western 
Australia, observing between 700MHz and 1.8 GHz, 
with an instantaneous bandwidth of up to 288 MHz. 
ASKAP is equipped with phased-array feeds (PAF; 
Hotan et al. 2014; McConnell et al. 2016), capable of 
simultaneously forming up to 36 independent beams, 
covering some 30 deg?. 

SPT2349—56, along with all 22 of the lensed SPT- 
SMGs in the ATCA program, were observed by the 
Rapid ASKAP Continuum Survey (RACS, McConnell 
et al. 2020), covering the sky south of --41deg dec- 
lination at a central frequency of 887.5 MHz, using 
903 individual pointings with 15-minute observations. 
The beam size at the location of SPT2349—56 is 
24" x 13". We retrieved the ASKAP image surrounding 
SPT2349—56 using the cutout server. At the declination 
of SPT2349—56 the achieved RMS sensitivity is 189 pJy. 
The RMS is similar in the ASKAP images around the 
other 22 lensed SPT-SMGs, although the actual sensi- 
tivity depends on proximity to other nearby bright ra- 
dio sources (see Appendix C). The SPT2349—56 ATCA- 
detected source is not cataloged in the RACS, but we 
find a 4.60 peak approximately 5" from the ATCA 
source (shown in Fig. 1). 


2.3. ALMA observations 


Extensive ALMA properties of SPT2349—56 sources 
B, C, and G have already been published (Miller et al. 
2018; Hill et al. 2020; Rotermund et al. 2021). Here we 
present several new ALMA observations (Table 2), sup- 
porting our measurements of line emission in the context 
of searching for AGN. 


1 The name means ‘shared skies and stars’ in the local indigenous 


language, Wajarri Yamatji. 
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ALMA Band-4 imaging (150 GHz) was obtained un- 
der three different programs in Cycles 3, 6, and 8, all 
targeting the brightest peak of the LABOCA source, 
and tuned to place CO(7-6) (Vrest = 806.652 GHz) and 
[C1](2-1) (Vrest = 809.34 GHz) in the upper sideband, 
and para-H30(211-202) (Vrest = 752.033 GHz) in the 
lower sideband. 

The Cycle 3 program 2015.1.01543.T (PI: K. Lacaille) 
was observed on March 20, 2016. The array was in the 
C36-2/3 configuration with baselines ranging from 15 to 
460m, and provided a naturally-weighted synthesized 
beam size of 0.88". Pallas and J2343—5626 were used 
to calibrate the amplitude and phase, respectively. The 
Cycle 6 program (2018.1.00058.S; PI: S. Chapman) ob- 
servations were obtained on 2018, October 3°¢ in the 
C43-6 array configuration with baseline lengths of 15 to 
2500 m, giving a naturally-weighted synthesized beam 
size of 0.28". J2056—4714 was used to calibrate the 
amplitude, while J2357—5311 was used to calibrate the 
phase. Lastly, the Cycle 8 program (2021.1.01313.S; PI: 
R. Canning) observations were obtained on 2022, July 
27. These observations used the C-6 array configura- 
tion with baselines of 15 to 2500 m, giving a naturally- 
weighted synthesized beam size of 0.27". J2357—5311 
was used to calibrate the amplitude, while J2336—5236 
was used to calibrate the phase. 

The Cycle 8 program (2021.1.01313.S) also observed 
CO(11-10) (Vrest = 1267.01 GHz) and continuum at 
about 230 GHz in Band 6. These observations, carried 
out on 2022, September 1, used the C-4 array configu- 
ration with baselines of 15 to 784 m, giving a naturally- 
weighted synthesized beam size of 0.47". J2357—5311 
and J2258—2758 were used to calibrate the amplitude, 
while J2357—5311 and J2336—5236 were used to cali- 
brate the phase. 

We also make use of previously-published Band 7 
(345 GHz) ALMA Cycle 5 and 6 observations (Hill 
et al. 2020). The deep 0.5"-resolution (ie. synthe- 
sized beam) Cycle 5 data contain the CO(16-15) line 
(Vrest = 1841.35 GHz) and an OH doublet; each of the 
doublets is actually composed of a triplet whose fre- 
quencies are about 0.01 GHz separated, which is com- 
pletely unresolved by our spectral resolution, so we 
consider the OH line to be a doublet. The mean 
frequencies of the doublet are Vrest = 1837.80 GHz and 
Vrest = 1834.74 GHz). These lines are present in the up- 
per sideband, which was not previously analyzed or pub- 
lished. The high-resolution Cycle 6 data described by 
(Hill et al. 2020) has a synthesized beam of about 0.2" 
and is here used to further analyze kinematics through 
a moment analysis of the [Cu] line (Section 3.3). 


All the data were calibrated using the standard 
observatory-supplied calibration script. Imaging was 
done using the CASA task tclean, using Briggs weighting 
with a robust parameter of 0.5, and in all cases channel 
widths were averaged down to a common 15.625 MHz. 
The Cycle 6 and 8 observations covering the CO(7-6), 
[C1] (2-1), and H20 lines were combined in uv space and 
then imaged together, while the Cycle 3 observation was 
imaged separately. We chose this approach as the two 
data sets did not overlap entirely in frequency, which led 
to artefacts in the imaging step. The higher-resolution 
Cycle 6 and 8 data cubes were then convolved to match 
the resolution of the Cycle 3 data (about 0.88"). The 
continuum was subtracted using the task imcontsub af- 
ter flagging all channels expected to contain line emis- 
sion based on previously-detected [C 11] lines given in Hill 
et al. (2020). At each spatial pixel, incontsub extracts 
a one-dimensional spectrum and calculates the average 
over all channels not flagged by the user, then subtracts 
this average and returns a continuum-subtracted data 
cube. 

The same apertures used by Hill et al. (2020) to ex- 
tract [C11] line strengths and 350-GHz continuum flux 
densities were applied to sources B, C, and G in order 
to extract one-dimensional spectra for each line. The 
Cycle 3 and Cycles 6+8 CO(7-6), (C 1|(2-1), and H20 
spectra were averaged to produce a final spectrum. De- 
tails on how line strengths and continuum flux densities 
(including our procedure for deblending lines) are given 
in Appendix B, and the spectra are shown in Figs. 9 — 
11. All new continuum flux densities and line strengths 
are listed in Table 3, and the new continuum measure- 
ments are also shown in Fig. 2. 


3. RESULTS 
3.1. Identifying and characterizing radio sources 


We first searched for radio sources at the positions 
of known ALMA and optically-identified members of 
the SPT2349—56 protocluster. There is one strong ra- 
dio detection at 2.2 GHz (5$955-—2144Jy) found near 
the SPT2349—56 core with ATCA (detected at 80), 
which corresponds to a less robust (4.60) detection 
with ASKAP at 888 MHz (Fig. 1 and Table 1). The 
ATCA source with a much smaller beam encompasses 
the bright central ALMA sources, named B, C, and G 
based on their rank-ordered 850 um flux densities Miller 
et al. (2018).? It is unclear from positional uncertainty 
and beam size whether the emission comes from all three 


2 These three sources are named C3, C6, and C13 in Hill et al. 


(2020) based on their rank-ordered [C 11] line strength. 
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galaxies or just a single source. Irrespective of this, 
the strong radio emission would be in excess from that 
expected from the far-infrared (FIR)-radio correlation 
(Helou et al. 1985). We analyse these issues in detail in 
Section 3.2. 

There are no other significant (>30) ATCA or 
ASKAP detections of any known protocluster members 
(Fig. 1). The FIR-radio correlation for star-forming 
galaxies (Helou et al. 1985; Ivison et al. 2010) would 
imply S2.2 2212 uJy for a Ss5o = 5 mJy source at z — 4.3. 
The ten brightest SPT2349—56 SMGs (excluding B, C, 
and G) span 0.8-15 mJy, with an average of 4.7 mJy. 
Thus even the brightest SMGs would only be expected 
to be at the 1c level in our ATCA map. A radio stack- 
ing analysis on these remaining ten brightest SMGs finds 
(11.0+10.0) uJy, which is completely consistent with the 
average 2.2 GHz emission expected from the FIR-radio 
correlation, ($55) =12yuJy. Stacking on all 40 known 
cluster members yields —5.0+5.8 uJ y. 

We then consider if there might be other radio sources 
in SPT2349—56 that could be cluster members. We 
searched for robustly-detected radio sources in the sur- 
roundings of SPT2349—56 out to 1 Mpc in projection 
(140" in radius) from the core, roughly the region stud- 
ied with ALMA by Hill et al. (2020). We find two ATCA 
sources above 5c, ID2 and ID3 in Table 1. ID2 is identi- 
fied to a bright star, and is also detected by ASKAP. 
ID3 has a clear optical counterpart, which does not 
have properties (especially non-detections in the g-band) 
of optical sources likely to be near z —4.3 (Rotermund 
et al. 2021). 

We thus focus on the properties of the central ID1 
radio source, starting with the positional uncertainty, 
Aa. From Condon (1997), we can derive the synthe- 
sized beam positional uncertainty for the ATCA and 
ASKAP detections, assuming that the beam is a single 
2D Gaussian with an RMS ‘width’ c =FWHM/2.354 
in each coordinate. In the limit where centroiding 
uncertainty dominates over systematic astrometry er- 
rors and for uncorrelated Gaussian noise, we have 
Aa =0.6 (SNR) FWHM. For both the ATCA and 
ASKAP sources in SPT2349—56, we have confirmed 
that the source size and position angle is indistinguish- 
able from other brighter, unresolved sources in the field, 
in agreement with the synthesized beam. We conclude 
that the SPT2349—56 radio source is unresolved with 
our current data. 

For the ATCA source (ID1) detected at SNR=7.9 and 
a beam size of 4" x 8" (PA — 27 deg east of north), the 
positional uncertainty is therefore 0.3" x 0.6". For the 
ASKAP source detected with SNR of 4.6 and a beam 
size of 24" x 13" (PA — 89 deg east of north) the po- 


sitional uncertainty is therefore 3.0" x 1.7". There is 
a 5.1” roughly northern offset between the ATCA and 
ASKAP sources, which is consistent at the joint 20 level. 
The ASKAP centroid is most consistent with the ALMA 
source A. Comparison of our wider field ATCA map 
and the ASKAP RACS map reveals that the major- 
ity of the sources show excellent astrometric alignment, 
but we also identify a few other ATCA sources with 
ASKAP counterparts with several arcsecond offsets (see 
Appendix A). In two cases, there is a robust association 
of the ATCA position to other cataloged objects (from 
2MASS), suggesting the offset to the ASKAP position 
is likely due to measurement error. For ID1, the more 
robust ATCA position and association to the B, C, and 
G galaxies in SPT2349—56 is the most likely interpreta- 
tion, with the ASKAP source being assumed to be en- 
tirely related to the AT'CA source for the purposes of de- 
riving a radio spectral index. The 5" offset is not entirely 
unexpected, but may be significant enough to require a 
physical interpretation rather than just measurement er- 
ror (Appendix A). It could for instance be related to a 
radio core-jet morphology. However, as noted, the 30m 
minimum baselines of ATCA would not resolve out flux 
on scales smaller than several arcmin. While it's not 
clear why the sources are offset, it appears more likely 
to be instrumental than physical based on the analysis 
in Appendix A. 


3.2. Physical interpretation of ID1 


We first constrain the radio spectral index to esti- 
mate and compare luminosities between sources. The 
radio source ID1 has a steep spectrum with an index of 
a= — 1.58 £ 0.28, constrained by the ASKAP 888 MHz 
detection, and the non-detections at 5 and 9 GHz. The 
uncertainty can be estimated by propagation of errors 
on the two frequencies as follows: 


V SNR;2 + SNRSÀ 
Aa = : (1) 
1n(2.2/0.89) 

In Appendix C, we describe a MCMC method to as- 
sess the uncertainty for any number of spectral measure- 
ments, and show this distribution in Figure 2. The spec- 
trum is too steep to be consistent with synchrotron ra- 
diation due to shock acceleration of cosmic ray electrons 
from supernovae (i.e., star formation), where Thomson 
et al. (2014) recently constrained a = —0.79+0.06 specif- 
ically for high-z SMGs. The steep SPT2349—56 spec- 
trum seems to demand an AGN interpretation. 

'The radio luminosity can then be assessed by assum- 
ing it is associated with one of the central SPT2349—56 
galaxies at z—4.3. With a specific luminosity of 
L3.3 = (4.4+0.3) x 10% W Hz- 1, it is far larger than ex- 
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Figure 3. Left: Redshift vs. 1.4 GHz radio power using the GOODS-N sample (Barger et al. 2017) and radio-excess candidates 
from the ALESS sample (Thomson et al. 2014). The detection threshold for the GOODS-N radio sample is shown with the blue 
curve. Red circles show sources detected above the 3c level at 850 um, while black circles show sources not detected at this level. 
SPT2349—56 has about 10 times more radio power than any radio source found in GOODS-N, although it has more than 500 
times lower radio power than other well-studied radio-loud galaxies that were used to identify high redshift protoclusters. Right: 
1.4 GHz luminosity over 350 GHz luminosity vs. redshift for the submm sources with spectroscopic redshifts in GOODS-N (red 
circles), and lower limits on radio sources undetected in the submm (red triangles). Also shown is the ATCA survey of lensed 
SPT-SMGs described in Appendix C (purple squares). The blue dashed line region shows where the submillimeter luminosities 
and radio luminosities produce consistent estimates of SFRs. None of the GOODS-N SMGs show any excess radio emission 
over the FIR-radio relation, while two of the ALESS sources do have a clear excess. Some of the higher redshift SPT-SMGs 
show a marginal radio excess, discussed in Appendix C. SPT2349—56 is about 100 times higher than the median relation at 
rest 1.4GHz. The HzRGs lie about 500 times higher with their measured Sagso = 6 — 12 mJy. 


pected from star formation through the FIR-radio cor- 
relation. For reference, the FIR-radio correlation for 
star-forming galaxies (Ivison et al. 2010) would im- 
ply Lo. —2.4 x 1074 W Hz-! for a similar Sg59 = 5 mJy 
source at z=4.3. Adopting the measured spectral 
index above, the radio excess increases to greater 
than a factor 100 at a rest-frame of 1.4GHz, with 
L1.4, rest = (2.40.3) x 1026 W Hz}. This strong radio 
excess suggests the presence of an AGN (e.g. Guidetti 
et al. 2017); however, the radio emission is still distinctly 
less luminous than powerful radio galaxies, like those re- 
siding in other structures studied at these redshifts, by 
a few orders of magnitude (Fig. 3). MRC 1138, for in- 
stance, is almost 1000 times more powerful in radio, and 
it is also hosted by the obvious BCG of the protocluster 
(e.g. Hatch et al. 2009). 

We then compare SPT2349—56 to radio sources from 
the literature. In Fig. 3, the redshift versus radio power 
is shown using the 0.3deg? GOODS-N VLA sample 
(Barger et al. 2017), which is highly complete in spec- 
troscopic redshift. We compute the rest-frame radio lu- 


minosity using the equation 


Lia = (4nd; 814/1079) (14 z) (ert) ergs ! Hz !, 
(2) 
where dz is the luminosity distance (in cm) and 
$14 is the flux density in units of uJy observed at 
1.4 GHz. This equation assumes 5$, x v^, and we adopt 
a radio spectral index of a= —0.8 (Ibar et al. 2010) 
for the GOODS-N sources, and the measured a for 
SPT2349—56 and the literature HzRG sources (in fact 
all very close to —1.6). Shown for comparison are 
several well-studied HzRGs that were used as beacons 
to uncover massive galaxy overdensities: MRC 1138 
(Large et al. 1981; Seymour et al. 2012); TN J1338 (De 
Breuck et al. 1999); and 8C 1435 (Lacy et al. 1994). 
SPT2349—56 has around 10 times more radio power 
than any radio source found in GOODS-N, but it has 

less than 500 times the radio power of these HzRGs. 
Figure 3 also directly assesses the departure of 
SPT2349—56 from the radio-FIR correlation by plot- 
ting the luminosity ratio of 1.4 GHz to 350 GHz versus 
redshift for all GOODS-N submm sources with spec- 
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Figure 4. Moment maps of B, C, and G. Left: 850-um continuum (red contours) from high-resolution (0.3) combined Cycle 


5 and Cycle 6 ALMA data (Hill et al. 2020), shown overlaid over 3-orbit HST F160W imaging (Hill et al. 2022). 


Mid-left: 


[C 11] moment-0 maps from Cycle 6 high-resolution ALMA data. Mid-right: [C11] moment-1 maps (in velocity units) from 
Cycle 6 ALMA data, with the zero velocity centered at the peak of each galaxy's [C 11] emission. All three sources show a 
clear velocity gradient (listed in Table 4, along with dynamical mass comparisons). Right: [C11] moment-2 maps (in velocity 
units) are shown from lower-resolution data to increase the SNR, revealing centrally concentrated dispersions. In all panels, the 
synthesized beam FWHM is shown in the bottom-left corner. 


troscopic redshifts (Barger et al. 2014; A. Huber in 
prep.). All of the submm sources in GOODS-N are 
radio-detected, even at z — 5.2, and the submm luminos- 
ity and radio luminosity produce consistent estimates 
of the SFRs for all sources — there is no sign of AGN 
from their radio emission. A similar analysis of the sub- 
set of gravitationally lensed SPT-SMGs also observed 
from ATCA in this program (Appendix C) suggests the 
majority (8796) also follow this relation; however, there 
are three very significant outliers in this sample which 
is most likely attributed to an AGN contribution from 
the foreground lensing galaxy (discussed further in Ap- 
pendix C). SPT2349—56 is an outlier by a factor of 
about 100 from this envelope (assuming the radio emis- 
sion is coming exclusively from ALMA source C). The 
HzRGs shown in the left panel of figure 3 have com- 
parable Saso = 6 — 12 mJy to other SMGs shown (e.g., 
Dannerbauer et al. 2014; De Breuck et al. 1999), and 
would remain about 500 times above SPT2349—56 in 
the radio/submm ratio plot in the right panel. By con- 
trast, the GOODS-N radio sources without submm de- 
tection rise significantly above this envelope, into the 
AGN regime. 

Thomson et al. (2014) have used deep JVLA (1.4 GHz) 
and GMRT (610 MHz) to study the 76 ALMA-identified 
SMGs in the CDFS field (the ALESS survey - e.g., 
Simpson et al. 2014). They find four SMGs whose radio- 
FIR values are > 20 above the sample median, which 
they classify as potential AGN. The most robust of these 
(ALESS 066.1) is a strong X-ray source with an inverted 
radio spectrum (a > 0.51). Of the remaining three, one 
(ALESS 014.1) has a flat radio spectrum (a > —0.1) 
and an obviously high radio luminosity, while the other 
two (ALESS 094.1 and ALESS 118.1) have spectral in- 
dex limits consistent with star formation (a ~ —0.8). 
We show these four SMGs in figure 3, where it is clear 


that none are comparable to SP T2349—56 in radio lumi- 
nosity or departure from the FIR-radio relation. In fact, 
two of the four are not at all unusual in their properties 
relative to the other samples. 

Radio emission provides an extinction-free probe of 
AGN (which even X-ray cannot claim, since practical 
sensitivity limits preclude the detection of the most ob- 
scured, Compton-thick AGN with Ng > 10?^cm-?). 
Traditionally radio AGN are divided into two sub- 
sets (Padovani 2017): (i) radio-loud AGN Li4 > 
1024 W Hz-!, which exhibit steep spectrum radio jets 
and lobes on kpc scales (Yun et al. 1999); and (îi) radio- 
quiet AGN, with flat-spectrum, lower luminosity radio 
emission, typically contained within a compact, several 
pc, core (Blundell & Kuncic 2007). SPT2349—56 is 
solidly a radio-loud AGN, whereas most of the other 
candidate AGN found in the surveys described above 
(GOODS-N and ALESS) cannot clearly be defined as 
such. 


3.3. Resolved properties of the ‘BCG’ sources 


Given the radio detection in SPT2349—56, it is of 
interest to assess the properties of the B, C, and G 
ALMA sources, and to compare them to other proto- 
cluster members. As noted, these are three very submm- 
luminous sources in the core region (Sgso =6.7 mJy for 
B, 4.7mJy for C, and 1.3mJy for G), with only source 
A being brighter, although two even more luminous 
sources are present in the northern extension (sources 
N1 and N2; Hill et al. 2020). 

The most distinguishing features of this trio (beyond 
their flux-ordered source names serendipitously spelling 
out ‘BCG’) are their locations near the center-of-mass 
of the cluster core, and their immediate environment. 
They are very close neighbours (they lie within an arc- 
second of each other), and are likely to be interacting. 
Further, there is a notable arc seen in [C 11] surrounding 


10 


CHAPMAN ET AL. 


Table 3. Continuum and line properties of B, C, and G. S147 and S231 are the continuum flux densities at 147 and 231 GHz, 
respectively, while the other columns provide various line strengths (the line is indicated by the subscript). The OH doublet 
arises from blended hyper-fine triplets centered at 1835 and 1838 GHz, and the H2O line is the para-211-2o2 line. 


ID RA,Dec S147 S231 Fcoüs-15) Fcoai-10)  FCo(r-6) Fazo Fon Ficy(2-1) 

ny ny Jykms”!  Jykms ! Jykms  Jykms ^ Jykms ^ Jykms! 
B  23:49:42.79, -56:38:24.0 5894-15 3322+156 0.12+0.05 0.26+0.09 0.73+0.05 0.22+0.02 0.93+0.11 0.46+0.03 
C  23:49:42.84, -56:38:25.1 336411 1810+118 0.10+0.06 0.1740.05 0.61+0.03 0.15+0.01 0.83+0.08 0.350.03 
G  23:49:42.74, -56:38:25.1 181423 1369 0.10+0.05 0.12+0.05 0.18+0.02 0.02+0.01 0.22+0.06 0.10+0.02 


Table 4. Physical properties of B, C, and G. SFRg,o is the SFR estimated using Eq. 5, Sgs0/Fu20 is the ratio of 850 um 
continuum flux density (from Hill et al. 2020) to H20 line strength, and SFRum/SFRan,o is the ratio of the FIR-derived SFR 
(from Hill et al. 2020) to the H>O-derived SFR. Vp-p is the peak-to-peak velocity from moment-1 maps, while FWHMcen is the 
central velocity dispersion (multiplied by 2v 21n 2) from moment-2 maps (see Fig. 4), and FWHMint is the width of the [Cu] 
line after fitting a single Gaussian to the lines shown in Fig. 9-11. Mayn, disk is the dynamical mass derived using Vp-p and a 
disk model (Eq. 3), while Mays, cen and Mayn, int are dynamical masses derived using the velocity dispersion measurements and 


Eq. 4. 
ID  SFRu,o Saso/FH;0  SFRyum/SFRgm,o0 Vp-p  Mayn,aisk FWHMcen Mayn,cen FWHMint Mayn, int 
Mo yr! 1072 km ^s kms! 10" Mg kms! 10? Mo kms ! 10? Mo 
B 11004410 3143 0.8493 600+50 18.2+2.2 5404-20 11.0+0.8 | 612410 14.0+0.5 
C 750+280 31:52 (gd 2404350 2.90.8 2804-20 2.9+0.4 35845 4.7+0.2 
G 80+60 652-23 agus 6904350  18.14+2.5 5204-20 7.6-+0.8 901454 22.8+42.8 


the three galaxies (Hill et al. 2020; N. Sulzanauer, in 
prep). Source C does distinguish itself with an anoma- 
lously narrow [C 11] and CO(4-3) line width for its lumi- 
nosity (Hill et al. 2020). Rotermund et al. (2021) iden- 
tified C as a significant outlier from the SPT2349—56 
galaxy sample in its Mayn/Mgas ratio inferred from the 
narrow CO(4-3) line width and large luminosity, similar 
to many high-z QSOs (e.g., Narayanan et al. 2008; Wal- 
ter et al. 2009; Hill et al. 2019), where selection effects 
favoring face-on orientation offer viable explanations. It 
is also noteworthy that source C has by far the largest 
stellar mass of any cluster member (> 10!! Mg Roter- 
mund et al. 2021; Hill et al. 2022). It is associated with a 
bright and very compact HST F160W source (Hill et al. 
2022), as shown in Figs. 1 and 4, and it has been sug- 
gested to be the seed of a growing BCG galaxy in this 
ongoing mega-merger (Rennehan et al. 2020). 


3.3.1. [Cu] kinematics 


We consider here a more detailed analysis of the kine- 
matic properties of the B, C, and G galaxies. Using high- 
resolution Cycle 6 [C 11] data (Hill et al. 2020), which has 
a synthesized beam of about 0.2”, we construct moment 
0, 1, and 2 maps of the B, C, and G sources and ana- 
lyze the resolved velocity and dispersion fields. We use 
the CASA task immoments, focusing on channels between 


+30 of the best-fit [Cu] line, and masking pixels < 4 
times the RMS per channel. Since second moments are 
particularly sensitive to noise (being a squared term), we 
use uv-combined Cycle 5 and 6 data cubes (described in 
Hill et al. 2020) to calculate the moment 2 maps; for 
reference, the resolution of the combined data is about 
0.3”. The results (moments 0, 1, and 2) are shown in 
Fig. 4. 

All three sources show a clear velocity gradient and re- 
solved, centrally-concentrated dispersion, characteristic 
of rotationally-supported disks. From these velocity gra- 
dients and velocity dispersion maps we extract peak-to- 
peak velocities, Vp-p, and central velocity dispersions, 
FWHMeen. We draw a line along the semi-major axis 
of each galaxy, then from the moment 1 map calculate 
the velocity difference between the two ends, and from 
the moment 2 map extract the velocity dispersion at the 
midpoint of the line. We find that moving the position 
angle of the line by +10 deg and moving the midpoint 
of the line by 5 pixels results in a peak-to-peak velocity 
change of 450 kms-! and a central velocity dispersion 
change of +10kms~! (+20kms~! in FWHM), so we 
quote these as our uncertainties. The results are given 
in Table 4, multiplied by a factor of 2V21n 2 to estimate 
a FWHM. 
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We use these peak-to-peak velocities and central dis- 
persions to estimate masses assuming a disk model, with 
the enclosed dynamical mass given by 


Mayn, dist[Mo] = 2.35 x 10° [V5 /(sin(i))] R, (3) 
where Vp-p is the peak-to-peak velocity in kms, 
R is the radius in kpc, and i is the inclination an- 
gle of the galaxy. We adopt a mean inclination suit- 
able for a collection of randomly oriented disks of 
(sin(i)) =a / 40.79 (see Law et al. 2009), and we use 
the half-light radii from Hill et al. (2022), estimated by 
fitting Sérsic profiles to the high-resolution ALMA [C 11] 
moment 0 images. The results are given in Table 4. 

The dynamical masses were derived previously (Roter- 
mund et al. 2021) from the unresolved velocity disper- 
sions, using the width of the integrated [C 11] lines shown 
in Figs. 9-10, with an assumption about the structure 
of the source based on the virial theorem, using the re- 
lation 


Mays [Ms] = 2.81 x 10° FWHM’ R, (4) 


where FWHM is a one-dimensional velocity dispersion 
(multiplied by a factor of 2/21n2) in kms~!, and R is 
the radius of the virialised structure. First, we use the 
resolved central velocity dispersion, FWHMeen, adopt- 
ing the [C 11] size measurements from Hill et al. (2022) 
and the central resolved velocity dispersions from the 
moment 2 maps (Table 4). Next we use the width of 
the integrated [C r1] line, FWHM;,;, obtained by fitting 
a single Gaussian model to the [C i1] spectra shown in 
Figs. 9-11 and given in Table 4, again using Eq. 4 and 
the same [CII] size measurements. The resulting dy- 
namical masses are provided in Table 4. 

Considered in the context of a disk model, source C 
does show a similar dynamical mass comparing both its 
central and integrated velocity dispersion (Table 4, and 
Rotermund et al. 2021); however, it still appears to have 
substantially lower mass (six times lower) than B from 
any kinematics analysis. Inclination is reasonably con- 
strained, since the aspect ratio of these galaxies is re- 
solved by ALMA. While it remains an uncertainty in 
any mass modelling, the aspect ratios of B and C are 
similar at ~1.8 (major to minor axis). 

Sources B and G have similarly large inferred disk 
masses (18 x 101! Mo). However, the distinct double- 
horned profile of source G (Appendix B) is direct evi- 
dence for a rotating disk or bar-like structure at high 
inclination (explaining the broad velocity profile), while 
the profile for B is possibly due to a tidal torque in re- 
sponse to the interaction with C. Source G also has a 
higher aspect ratio (2.3, major/minor axes) in moment-1 


than B and C, suggesting the disk is seen closer to edge- 
on. Explicitly using this higher implied inclination in 
Eq. 3 brings down the disk mass estimate by 2596, more 
consistent with the much lower gas mass of G compared 
with B and C. 

It is noteworthy that in projection at least, B is 
counter-rotating relative to C. Several studies have pre- 
dicted that mergers configured with counter-rotating gas 
disks should lead to the most intense starbursts, and 
conditions for fueling the SMBHs (e.g., Mihos & Hern- 
quist 1994, 1996; Di Matteo et al. 2007; Salomé et al. 
2012). 


3.3.2. Submm line properties 


We then consider line diagnostics to elucidate which 
of the three might be most likely to host the radio-AGN. 
We first assess the [Cu]/FIR ratio, which has been 
shown to highlight AGN with a deficit compared with 
star-forming galaxies (e.g. Stacey et al. 2010). However 
at high luminosities, both AGN and SMGs (without ob- 
vious AGN) exhibit similar deficits in the ratio. Hill 
et al. (2020) have shown that all three of B, C, and G 
are ‘deficit sources’ in [C 11]/FIR, inhabiting similar re- 
gions in the [C r]/FIR-to-FIR plot as many luminous 
AGN. However, this work also showed that all 12 of the 
most luminous SMGs in SPT2349—56 have comparable 
(C u]/FIR ratios, and none of these are obviously AGN 
from any available diagnostics. 

One possibility to consider is that the FIR estimates 
are being affected by an AGN in one of B, C, or G. 
Since the shortest wavelength measured by ALMA is 
160um in the rest frame, the peak of the SED is not 
sampled, and there is little constraint on whether the 
dust might be substantially hotter than the Tg ~40K 
estimated in Hill et al. (2020). To test this we make use 
of the para-H2O0(211-202) lines observed in the ALMA 
Band 4 dataset (Table 3 and Figs. 9-11). H2O is 
strongly coupled to the FIR radiation field whether it 
is being produced by star-formation or AGN (Omont 
et al. 2013). Jarugula et al. (2021) compiled a sample 
of low- and high-z submm galaxies with para-H5O0(21,—- 
202) measurements (including two sources, SPT0346—52 
and SPT0311—58, from the same parent sample as 
SPT2349—56), and found that a simple single-parameter 
scaling relation described the correlation between Lq,o 
and SFR (derived from FIR) of the form 


SFR [Mo yr 1] = (2.07 + 0.75) x 10 ?Lu,o [Lg]. (5) 


A simple test is to first take the ratio of 850jum contin- 
uum flux density to H2O line strength, where measure- 
ment errors are mostly small. These values are given in 
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Figure 5. SLED for B, C, and G in the CO Jupper = 2, 4, 7, 11, and 16 transitions (new data for the Jupper = 7, 11, 16 lines are 
shown in Appendix B). The ATCA detection of CO(2-1) (Miller et al. 2018) is shown as an upper limit as it is an unresolved 
measurement over the core region including at least B, C, and G. Open circles illustrate the J = 2 division if the luminosities 
scale from J = 4. We compare to the Lrrr —3 x 10? Lo AGN-dominated galaxy Mrk231 (van der Werf et al. 2010), and the 
Ly —3 x 101 Lg starburst M82 (Kamenetzky et al. 2012), here normalized to Mrk231 at CO(7-6). In the right panel the 
SPT2349—56 galaxies are also normalized to the CO(7-6) luminosity of Mrk231 for comparison of the excitation curves. Source 
G is undetected in the Jupper = 11 and 16 lines, while B and C (shown as 20 upper limits) are marginally detected in CO(16-15). 


Table 4, where we have used Sg50 values from Hill et al. 
(2020). Using the same modified blackbody SED as in 
Hill et al. (2020) to model the continuum flux density 
emission, a dust temperature of 40 K at a redshift of 4.3 
means that Sg59 =1 mJy corresponds to 115 Mọ yr7!, 
and so Eq. 5 implies Sg59/ Fu,o = (42 x 1073) km ^! s; B 
and C sit significantly below this value, implying they 
might have higher FIR than currently estimated. In con- 
trast, G is significantly above the relation, which could 
imply cooler dust and lower FIR than previously esti- 
mated. This would make G less of a deficit source in 
[Cu] /FIR, and less likely to be considered an AGN by 
this criterion. 

Using Eq. 5, we can also estimate SFRs directly for 
B, C, and G using our measured para-H2O(211—202) line 
strengths (Table 4), and compare these with the SFRs 
from the FIR (taken from Hill et al. 2020). These re- 
lations show the same behaviour as our ratio of mea- 
surements above, that B and C both may have higher 
FIR (from hotter dust) than estimated from our current 
ALMA data. We note that systematic errors in convert- 
ing to these physical quantities are large (as listed in 

We next consider the CO spectral line energy distribu- 
tion (SLED), which can distinguish AGN with high ex- 


citation lines driven by X-ray dominated regions (XDRs, 
e.g. van der Werf et al. 2010). Here we present higher- 
J CO transitions (Section 2) than have previously been 
published in Miller et al. (2018). All of B, C, and G are 
well-detected in CO(7-6) observations, while B and C 
are detected in CO(11-10). Remarkably B and C may 
be marginally detected in CO(16-15) due to the sensi- 
tivity of the deep Band-7 ALMA data presented in Hill 
et al. (2020), although strictly they are upper limits. All 
the new line channel maps and one-dimensional spectra 
are shown in Appendix B. SLEDs are shown in Fig. 5 
from the available J —2, 4, 7, 11, and 16 transitions, and 
compared to AGN and starburst templates. The ATCA 
detection of CO(2--1) (Miller et al. 2018) is shown as 
an upper limit as it is an unresolved measurement over 
the core region including at least B, C, and G. However 
for illustration, we show the division into B, C, and G 
of the integrated J — 2 luminosity assuming they scale 
with the J — 4 fluxes. None of B, C, or G appear to have 
high excitation SLEDs similar to AGN like Mrk231 (e.g. 
van der Werf et al. 2010), and are all similar to or less 
excited than M82 at high-J (Kamenetzky et al. 2012), 
with the caveat that G has only upper limits beyond 
J= 
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Source B has the highest excitation SLED confirmed 
of the three, lying near the M82 SLED. Source B also has 
a stronger cool/warm gas component similar to Mrk231. 
However, all three are reasonably characterized with a 
combination of cool and warm star-forming Photo Dis- 
sociation Regions (PDR) components, and without sig- 
nificant XDR contributions. Detailed SLED modelling 
of SPT2349—56 sources will appear in a future contri- 
bution. 

Finally, the 1634m OH doublet in B, C, and G can be 
compared to Runco et al. (2020), who studied 178 local 
galaxies in six of the 14 OH transitions in the FIR range. 
They found the highest frequency OHi63um (detected in 
25 galaxies) is the only OH doublet which is always in 
emission, with most transitions often appearing in ab- 
sorption. Runco et al. (2020) presented the correlations 
of the equivalent width, EW(OH), with various galaxy 
properties and line ratios, finding EW(163um) is not 
well established as a direct AGN indicator. For example, 
while galaxies with lower X-ray luminosities exclusively 
have low EW(OH), the full range of EW is seen for the 
highest X-ray luminosities (Runco et al. 2020). How- 
ever, a strong correlation is found for EW(OH) with the 
ratio of AGN activity to SFR, suggesting this is a bet- 
ter predictor of EW(OH) than the total AGN power. In 
figure 6, we compare the equivalent width, EW (163um), 
to local starbursts, LINERs and Seyfert galaxies from 
Runco et al. (2020). The EW(163um)=0.076um and 
0.095um measured for B and C respectively are amongst 
the highest found locally. (G is similarly high, but is 
only marginally detected in OH). Their EW(163jim) 
are more similar to values in local Seyfert galaxies than 
starburst galaxies, the latter having EW(163um)~0.02- 
0.05um. It is not yet clear at z > 4 what is “normal” 
for EW (163um), since the line has only before been de- 
tected locally. However, these results may provide ini- 
tial evidence that B and C do in fact present as AGN 
through some submm-wave diagnostics. 


3.3.3. Optical and near-infrared properties 


Finally, we summarize optical and near-infrared spec- 
tra taken with the Gemini and VLT observatories. A 
VLT XSHOOTER spectrum (Aobs = 0.35-2.4 um) was 
obtained which targeted B and C (Rotermund et al. 
2021), covering redshifted Lya through [OII]3727._ No 
lines were detected. A VLI-MUSE spectral cube was 
used to extract one-dimensional spectra at the locations 
of each of B, C, and G (Y. Apostolovski in prep.), but 
no lines are detected in any of these galaxies (nor any of 
the SPT2349—56 SMGs). 

Non-detections are not particularly surprising given 
the faintness of the galaxies at 6450 A, the wavelength 
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Figure 6. The equivalent width of the OHi¢63,m doublet 
versus IR luminosity. B, C, and G (blue circles) have similar 
values, although the error in G is large. We compare to 
all local galaxies detected in OHi63um from the compilation 
in Runco et al. (2020). Galaxies are colour coded by their 
classification as starbursts (lime), LINERs (green), Seyfert-1 
(red), and Seyfert-2 and intermediate types (orange). B, C, 
and G all appear well above local starbursts, although all 
three AGN types have a few examples this high. 


of redshifted Lya, where C has So.ssum = 0.061 uJy, 
while B and G are undetected to So63um < 0.01 pJy 
(Rotermund et al. 2021; Hill et al. 2022). The diffi- 
culty of spectroscopy in the near-infrared at 19,768 A 
in the vicinity of the redshifted [OII]3727 means these 
limits on line equivalent widths are also not particularly 
constraining. Nonetheless, strong AGN often exhibit de- 
tectable high excitation lines in optically faint, obscured 
SMG hosts (Chapman et al. 2003; Chapman et al. 2004; 
Chapman et al. 2005; Danielson et al. 2017), and it is 
surprising that this AGN in SPT2349—56 eludes all op- 
tical and near-infrared spectral detection. 


3.3.4. Concluding remarks 


Thus while the radio observations do not have suffi- 
cient spatial resolution to uniquely identify one of the 
three galaxies as the AGN, the source properties them- 
selves suggest source C could be a likely host, consid- 
ering mainly its large stellar mass, along with narrow 
emission lines, and high EW(OH). Sources with sim- 
ilar radio luminosities in the local Universe are typi- 
cally found in massive hosts. However, the more FIR- 
luminous and much more dust-obscured source B might 
also be a possible host for the AGN, given the large dy- 
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namical mass from kinematic modelling and the higher 
excitation SLED. Of course all three SMGs could have 
AGN components at the same time. 'The fact that they 
are likely strongly interacting dispels the typical duty- 
cycle arguments that would disfavor this scenario. To 
make progress, we will need deeper radio data with 
better resolution, and sensitive infrared spectroscopic 
observations now possible with the James Webb Space 
Telescope. 


4. DISCUSSION 
4.1. Inferring AGN properties from radio power 
4.1.1. Jet power and energy input to the ICM 


Radio jets are thought to provide an important feed- 
back mode in galaxy clusters by preventing the cooling 
of hot (X-ray) gas surrounding central galaxies (e.g., Mc- 
Namara & Nulsen (2012)). This is named “jet-mode” 
feedback and is associated to radio sources characterised 
by radiatively-inefficient accretion. However, radio jets 
can also drive massive gas outflows on galactic scales, 
another signature of AGN feedback. 

A theoretical relation between radio luminosity and 
radio jet power was determined by Willott et al. (1999), 
and can be used to estimate the kinetic energy output of 
AGN (e.g., Hardcastle et al. 2007). The jet power can be 
estimated by assuming that the mechanical power of the 
jet can be approximated as the energy of the detected 
radio cavity averaged over some timescale (e.g., Birzan 
et al. 2004). The X-ray-detectable “cavities” that result 
from AGN jet activity (O’Sullivan et al. 2011) allow us 
to quantify the heating experienced by the intra-cluster 
medium (ICM). The energy contained in these cavities 
comes from the product of the pressure and volume (pV) 
over the cavity. This is the work done by the jet to cre- 
ate the cavity, and the internal energy of the radio lobes. 
Under the assumption that the cavity is dominated by 
relativistic plasma, this becomes 4pV. Dividing the en- 
ergy of the cavity by the cavity age gives the power, 
Poay- 

Thus the most direct inference we can make from 
the radio properties of SPT2349—56 adopts a relatively 
tight correlation observed between radio power and cav- 
ity power (Cavagnolo et al. 2010; O’Sullivan et al. 2011; 
Panessa et al. 2015), where a fitted relation follows: 


log Peay = (0.35 + 0.07) log Li.4 + (1.85 +0.10), (6) 


yielding Pray = (3.3 + 0.7) x 10% W. This is strictly a 
lower limit to the jet power, and therefore energy injec- 
tion into the ICM. The true jet power depends on how 
the radio cavity is inflated (as described in Nusser et al. 
2006), with some energy from the jet being carried away 


by shocks. The relation of Li 4 to Peay is still affected 
by uncertainties due to the assumption that the cavity is 
dominated by relativistic plasma and the detectability 
of cavities within the sample used in O’Sullivan et al. 
(2011), as discussed in their work. 

This jet power is a sizeable amount of energy, given 
the potential well of the 101% Mg SPT2349—56 halo 
(see below) constrained from the central velocity disper- 
sion and radial distribution of cluster members (Miller 
et al. 2018; Hill et al. 2020). This is also a significant 
addition to the already abundant energy injection from 
the 6600 Mg yr! of star formation being experienced 
by the core of SPT2349—56 from summing the SFRs of 
all member galaxies found in these works (Miller et al. 
2018; Hill et al. 2020; Rotermund et al. 2021). We take 
the instantaneous injection of energy at z = 4.3 as 


Exin = 5 Mow”, (7) 
where Mout is the total amount of gas ejected per unit 
time by galaxies and v is the outflow velocity. While v 
is not measured in SPT2349—56 galaxies, the average 
outflow in high-z SMGs and other starforming galax- 
ies has been constrained with increasingly large samples 
(e.g., Banerji et al. 2011; Forster Schreiber et al. 2014). 
We adopt a typical 500 km s^! wind speed for the SN- 
driven outflows in each SPT2349—56 galaxy. A mass 
outflow rate can then be found by converting SFRs into 
mass outflow rates Mout by multiplying by a conser- 
vative mass loading factor n= Mout /SFR=1. m could 
even be greater than one based on observational (i.e., 
Newman et al. 2012) and theoretical work (i.e., Hop- 
kins et al. 2012). However, the same amount of metals 
is found in stars and the ICM, which suggests equal- 
ity, Mout SFR, (e.g., Renzini & Andreon 2014). We 
therefore obtain Euin = (3.3 + 1.1) x 10% W, where the 
uncertainty reflects both the range in SFR estimates and 
the range of likely wind velocities. This energy injection 
is remarkably similar to that found from the radio-loud 
AGN above from equation 6. 

Estimating the total mechanical energy injected by 
the radio jets requires an estimate of the radio source 
lifetime. Brienza et al. (2017) and Hardcastle et al. 
(2019) have suggested that remnant sources fade rapidly, 
with most of the observed remnant radio galaxies being 
relatively young, with ages between 50 to 100 Myr. 

With 7—100Myr for SPT2349—56, we find 
Emech = Peay X T= (1.0 + 0.2) x 10°4J, assuming only 
the uncertainty in the Peay scaling relation. 


4.1.2. Binding energy of the halo gas 


We then turn to estimating the binding energy of the 
gas in the SPT2349—56 halo. Giodini et al. (2010) 


A RADIO-LOUD AGN IN SPT2349 15 


demonstrated that the mechanical energy from jets is 
comparable to the binding energy (Epinding) in galaxy 
groups, while it is lower by a factor of 102-105 in clus- 
ters. Since the SPT2349—56 halo mass is comparable 
to a large group today, and the entire protocluster is 
expected to form a massive cluster by z = 0, it is thus 
of interest to investigate how Finding compares to our 
estimate of Emech- 

We define the binding energy as the total potential 
energy needed to push the ICM gas within Rsoo (the 
radius where the mean dark matter halo density drops 
to 500 times the critical density) beyond Raoo (the radius 
where the mean dark matter halo density drops to 200 
times the critical density, which we assume to be equal to 
the virial radius). Hill et al. (2022) estimated M200 (the 
mass contained within R99) to be (9+5) x 10? Mo, 
corresponding to A599 = (120+ 70) kpc at z —4.3, and 
Hoo can be computed if one assumes a density profile 
for the dark matter. 

Following Giodini et al. (2010), the binding energy is 
computed as 


Megas,500 
Evinding =f le(r) $(R200)] d Mgas 
R500 
=4 rf P(r) pzas(7) r° dr, (8) 
0 


where the constant term (R200) is small compared to 
the potential within R5oo and can be ignored, and pgas 
is the gas mass density. 

Assuming the gas mass density follows the dark matter 
density but scaled by a single gas-mass fraction param- 
eter, feas, we can adopt an NFW dark matter profile to 
write the binding energy as (see Giodini et al. 2010 for 
details) 


C500 l 1 
Epinaing = feast T Perit bc A r3 | a ER 2) dx (9) 


o (tex)? C 


where z — r/r,, with r, being the characteristic ra- 
dius related to the halo concentration parameter by 
c= R200/rs; c is a numerical factor that depends only 
on the halo concentration parameter c, A scales with 
Ms»ogo and also depends on c, and perit is the criti- 
cal density at the redshift of interest (here 4.3). We 
compute the concentration parameter using the mass- 
dependent relation of Macció et al. (2007); they find a 
linear trend between log c and log My, (which we assume 
is equal to Mago), and we find c— 7.9, corresponding 
to rs = 15kpc. We note that c— 5 is typically adopted 
for massive clusters > 10!^ Ma. With the concentra- 
tion parameter known, we calculate R500 = 80kpc and 
C500 = Rso0/Ts = 5.3. 


We cannot estimate the halo gas mass directly in 
SPT2349—56, beyond summing the measured cold gas 
masses in individual galaxies from the core region and 
inferring additional cool and warm gas components in 
the halo. Summing the Hə gas masses from the 23 
SMGs within the cluster core from Hill et al. (2020) 
yields Mgas,coo! = 3 x 1011 Mg. The unseen gas compo- 
nents in the halo are more uncertain. A trend observed 
in groups and clusters is an increase of the fraction of 
hot gas with total system mass (Connor et al. 2014), 
approximately following fgas x M°!~°-?, where 101? Mg 
groups typically have feas of around 10%. 

The Lx-M relation has been shown to remain approx- 
imately self-similar out to z 2 2 (Mantz et al. 2018), in- 
cluding X-ray-detected clusters at z —2 (Gobat et al. 
2011). However, for low-mass systems the gas mass frac- 
tions may evolve with redshift (Connor et al. 2014). Re- 
gardless, this in itself does not constrain the ICM gas 
fraction, which requires more detailed X-ray properties 
than Lx to be detected. 

Based on the above, we will assume for SPT2349—56 
a gas mass of 10% of the halo mass, or 9 x 10!! Mo, 
which nominally requires that Maas noc =6 x 1011! Mo, 
unless there are substantial cold flows feeding the 
submm galaxies (Dekel et al. 2009). We estimate 
Ebpinding = (1.5791) x 1054 J, where the uncertainty has 
been propagated from M»oo and the uncertainty in the 
Mooo-c scaling relation using a Markov chain Monte 
Carlo (MCMC) approach. 

The radio feedback alone therefore conceivably pro- 
vides all of the energy required to unbind the total gas 
in the cluster core. The stellar feedback has a compara- 
ble energy input, and could also be unbinding the cluster 
gas. However, the total energy is a minimum condition; 
the energy must also couple efficiently to the ICM. An 
energetic jet may not couple to the bulk of the ICM gas 
(Babul et al. 2013; Yang & Reynolds 2016; Cielo et al. 
2018). 

Any hot ICM established at z>4 may not be in hy- 
drostatic equilibrium since cold inflows likely dominate 
the flow of gas in protocluster halos (Dekel et al. 2009). 
The infalling gas only increases the energy required to 
inflate a bubble in the nascent ICM, acting as an addi- 
tive term to Ebinding. While L4 is fixed, the work done 
on an inflowing medium will be higher than for an ambi- 
ent static medium. Therefore the Peay (c 4pV) to Lua 
relationship might not hold when inflows dominate the 
halo. Yajima et al. (2022) and Trebitsch et al. (2021) 
have begun to explore some of these issues in hydrody- 
namical simulations of protoclusters, aiming to better 
understand AGN feedback and the impact of massive 
starburst galaxies in forming clusters. We leave more 
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detailed calculations to future work (D. Rennehan, in 
prep.). 
4.1.3. Inferred X-ray luminosity and accretion rate 


A correlation also exists between radio power and X- 
ray luminosity (Lx) for radio-loud AGN (Ballo et al. 
2012), although there is substantial scatter in this rela- 
tion. While the correlation appears to be similar over 
a large range (nine orders of magnitude) in X-ray lumi- 
nosity, there is a range of over 100 in Lx for a given 
radio luminosity in well populated areas of the correla- 
tion. The relation plotted in Ballo et al. (2012) is char- 
acterized at 5 GHz rest frame, which we measure almost 
directly (through the ASKAP detection). Using the cor- 
relation, we find that L; =7 x 10% W Hz-! in the radio 
corresponds to Lx = 1038 W (where the X-ray luminos- 
ity is between 2 and 10 keV). We conclude that the X-ray 
emission from the central AGN in SPT2349—56 can be 
easily detected by XMM- Newton or Chandra under the 
full range of possible Lx = 1037-39 W suggested by this 
correlation. 

Finally, taking source C as the most likely host, we 
can infer the SMBH mass from the stellar mass that has 
been well characterized for C (Rotermund et al. 2021; 
Hill et al. 2022). For M* 24 x 1011 Mo, the SMBH mass 
is 7 x 108 Mg (e.g. Ding et al. 2020). From this, we can 
infer the range of Eddington luminosities with respect to 
the range in X-ray luminosity constrained by the radio 
power. In other words, how close to the maximal rate of 
accretion is the SPT2349—56 AGN if its SMBH is close 
to that implied by the stellar mass of source C. In partic- 
ular, following Ballo et al. (2012), our measurements of 
Ls auz/Mpg constrain Lx / Lgaa to the range of roughly 
0.005 to 0.05, based on their distribution shown (their 
figure 10). Directly measuring the X-ray properties of 
SPT2349—56 will allow substantial progress in charac- 
terizing the system and its environment. 


4.2. Implications of the steep spectrum 


For an optically-thin synchrotron source, the spec- 
trum will steepen in spectral index from low to high 
frequencies by Aa — —0.5 if the source lifetime is greater 
than the timescale for energy-loss from the radiating 
electrons. This leads to a concave spectral shape with 
a characteristic bend frequency, v» (Kellermann et al. 
1969). Thus the age of the electron population within 
radio jets contributes to the steepness of the spectrum. 
Three effects will then decrease m, as the source red- 
shift increases (Krolik & Chen 1991): (1) for a fixed 
bend frequency v+» in the rest frame, the observed bend 
V, = v, / (1+ z); (2) losses due to inverse Compton scat- 
tering off the microwave background rise with redshift 
as (1 + z), so that for a fixed time electrons spend in 


the radiating region, the lowest energy electron that can 
cool has a frequency (or energy), which decreases with 
increasing redshift; and (3) flux-limited samples result 
in a selection effect that favors low v, at high-z. Sources 
must have higher emissivity at higher redshift to be in- 
cluded in the sample. They also must have stronger 
implied magnetic fields, and therefore more rapid syn- 
chrotron losses. 

A combination of these effects has been used to explain 
the observed trend that higher-redshift radio galaxies 
have steeper spectral indices (Carilli & Walter 2013; van 
Breukelen et al. 2009). The ultra-steep spectral indices 
of HzRGs (up to the a = 1.6 we find in SPT2349—56) is 
a main selection criterion for identifying these powerful 
radio sources in the distant Universe (De Breuck et al. 
2000; Broderick et al. 2007). All three HzRGs shown 
in Fig. 2 in fact have o very close to 1.6. It is of note 
that SPT2349—56 would have been discovered by these 
HzRG surveys over one to two decades ago had the radio 
source been 10-100 times more radio luminous, and even 
cursory submm followup would then have revealed the 
extended Sgs59 = 110 mJy source that belies its nature 
as a submm-luminous protocluster. 

A steep spectrum generally argues for self-absorbed 
synchrotron, and a lack of electron injection (e.g., Rad- 
cliffe et al. 2021). Thus the steep a in SPT2349—56 
could represent a dying radio source. In this case, the 
ATCA flux should be extended over the same area as 
the ASKAP data. Thus if there is no physical offset 
between ATCA and ASKAP, and this is just a mea- 
surement uncertainty, SPT2349—56 could be a young 
and completely unresolved compact radio source. On 
the other hand, this might be a “contained” or “frus- 
trated” radio source inside a dense medium, sometimes 
referred to as a compact steep spectrum source, or CSS 
(Padovani 2017), but an issue with this interpretation 
is that the luminosity of the source is low relative to 
these typical GHz-peaked sources. If self-absorbed syn- 
chrotron is contributing to the steep spectrum, the ob- 
servational constraints would mean that the break fre- 
quency is well below about 5 GHz in the rest frame. In 
principle this break frequency can provide a constraint 
on the age of the radio source, but since we do not con- 
strain this break with the current data, we do not pursue 
this further here. However if the radio emission is due 
to a CSS then it would have to be older than 500 Myr 
to have a break frequency below 0.9 GHz (4.8 GHz rest) 
(Padovani 2017). 


4.3. Connection to the LAB 


The powering sources of Ly-a blobs (LABs) have of- 
ten been identified broadly with the photoionizing emis- 
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Figure 7. Does the radio-loud AGN power the LAB? The 
background shows HST F160W imaging with ATCA 2.2 GHz 
contours (cyan) and the Y. Apostolovski et al. (in prep.) 
MUSE Lya contours (lime). The center of the LAB lies 
4.5" (31 kpc) from the radio source centroid. ALMA 850um 
contours are shown (coral), but sources M and N are too 
weak to see in this representation. The LAB is centered 
near the SMG, N, which was originally identified through 
its [C 11] emission and undetected in continuum (Miller et al. 
2018). N has Simm = 0.18mJy and an Lig = 4 x 10!! Lo, 
with an implied SFR=35 Mo yr |. 


sion from a close ionizing source (e.g., a QSO, Geach 
et al. 2009; Overzier et al. 2013), shocks (e.g., Taniguchi 
& Shioya 2000), or “cooling radiation” during gravita- 
tional collapse of the gas (e.g., Haiman et al. 2000). The 
SPT2349—56 LAB (shown in figure 7) was originally 
hypothesized to be heated by some combination of the 
three ALMA sources that reside near or within it (Y. 
Apostolovski et al. in prep.). However, given that the 
LAB center is only 4.5" (31 kpc in projection) offset from 
ALMA source C, it could instead be heated by the radio- 
loud AGN. The LAB is centered on the weak SMG, N, 
which was originally identified through its [C r1] emis- 
sion (Miller et al. 2018). N is a luminous infrared galaxy 
(LIRG) with Sasso =0.27+0.04 mJy, Lei =4 x 10H Lo, 
and a substantial M* —3x 101? MS. It is a plausi- 
ble, but somewhat unlikely power source for the lumi- 
nous LAB (whose total luminosity is 3 x 10% ergs-!, or 
3 x 10% W); source N fails to provide the necessary UV 
ionizing photons by at least a factor of ten, scaling from 
its meager R-band flux density of 0.37 uJy (similar to 
the analysis in Y. Apostolovski et al. in prep.). We can 
directly estimate the AGN X-ray emission expected for 
powering the Ly-a blob following Overzier et al. (2013), 


assuming that the fraction of ionizing photons that will 
cascade to Lya is 68% (case B recombination). This 
number likely exceeds the actual amount of ionizing ra- 
diation available due to the absorption by dust by a fac- 
tor of around 10, which we account for here. We then 
assume a radio-quiet QSO spectrum given by Richards 
et al. (2006). The predicted observed frame (0.2-12 keV) 
X-ray luminosity would be 2 x 1057 W, which is compa- 
rable to the low end of the expected range of Lx from 
the SPT2349—56 radio source, as discussed above. The 
radio AGN may therefore be at least as plausible a heat- 
ing source as N. 

Regarding the Lya blob being spatially offset from 
the AGN position, we note that in the radio source 
B3 J2330 at z 23.1 (Matsuda et al. 2009), the peak of 
the Lya emission was also found to be similarly offset 
from the HzRG itself. Even in the z=4 Distant Red 
Core (DRC) LAB, there is a roughly 3" (21 kpc) offset 
from the X-ray-emitting AGN that is proposed as the 
LAB's power source (Vito et al. 2020). However, these 
are rare cases. Venemans et al. (2007) showed that gen- 
erally the AGN is very near the center of the Lya halo, 
which grants some geometrical credence to the idea that 
the Lya halo is ionized by the central AGN's photons. 
In SPT2349—56, this is harder to argue, but the Lya 
could be completely absorbed by the copious amounts 
of dust in the core. The SPT pre-selection (as with the 
Herschel selection of the DRC) may favor finding sources 
with such offsets. 


4.4. AGN fractions in protoclusters 


As described in Section 3.1, with 27uJy RMS 
at 2.2GHz, we are sensitive to moderately-luminous 
and heavily-obscured z= 4.3 AGN among the 30 sub- 
millimeter galaxies identified in the SPT2349—56 struc- 
ture. They need to lie approximately 5 times above 
the radio-FIR relation to be significantly (50) detected 
by ATCA. In GOODS-N (Fig. 3), there are seven radio 
sources (all lacking submm detection) that satisfy this 
threshold, all of which lie at at redshifts less than 2. 
Another seven such radio sources lie 2-3.5 times above 
the relation, extending to a redshift of about 4, which 
would not be detected by our observations. The fact 
that all submm-detected sources in GOODS-N, and 74 
of 76 SMGs in ALESS, are consistent with the radio- 
FIR relation does signify that radio-loud AGN are not 
common amongst the submm-luminous population. No 
significant radio emission is found from any other (non- 
SMG) cluster members or candidates. With our cur- 
rent radio depth, the radio-AGN content among SMGs 
in this protocluster is constrained to be less than 1096 
(three of 30 members), and most likely 396 (assuming 
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C is the host of the ATCA radio source). However, the 
radio-loud AGN are only about 10% of the total AGN 
population in the field (Barger et al. 2007; Radcliffe et al. 
2021). The X-ray AGN fraction remains unconstrained, 
and given that many X-ray AGN are not radio emit- 
ters (Barger et al. 2007), our AGN fraction estimates in 
SPT2349—56 are lower limits. 

In the z 2 4.0 DRC protocluster (Oteo et al. 2018), 
a central galaxy is radio-undetected, but is a Compton- 
thick X-ray AGN. Only one of the three X-ray-identified 
AGN is detected in the radio - DRC6 (S5.5 = 128 wJy, 
$9 — 120 uJy), indicating a flat-spectrum source. In this 
case, the radio-AGN in the DRC lies towards the edge of 
the projected distribution of SMGs (offset from the core 
of the cluster). Thus without X-ray data, we cannot 
tell if the total AGN fraction of SPT2349—56 is differ- 
ent from that in the DRC (2396, Vito et al. 2020). As 
another example, in the core of the z = 3.09 SSA22 pro- 
tocluster, the SMGs have a 5096 X-ray AGN fraction, 
with four of eight SMGs detected by Chandra (Umehata 
et al. 2019), significantly larger than the DRC. 


4.5. Radio sources and cluster evolution 


Given that SPT2349—56 is conceivably the most mas- 
sive and active halo we know of at z » 4, an open ques- 
tion concerns the feedback or radio mode that this AGN 
is operating in, and how it is shaping the early core evo- 
lution of the cluster. With the current data, having 
only the two photometric points characterizing the ra- 
dio emission, and not even localizing it uniquely to one 
galaxy, we cannot definitively address these issues. Most 
radio-loud AGN appear to be hosted in recent or ongo- 
ing mergers (e.g., Ramos Almeida et al. 2012; Chiaberge 
et al. 2015). In this light it may not be too surprising 
to find a radio-loud AGN in the core of SPT2349—56. 
Given that the radio luminosity of SPT2349—56 is mod- 
est for an HzRG, we may be seeing a radio-loud AGN 
fueled via radiatively-inefficient flows with low accretion 
rates (Best & Heckman 2012). In this picture, the gas 
supplying the radio galaxy is frequently associated with 
hot X-ray halos surrounding massive galaxies, groups 
and clusters, as part of a radio-AGN feedback loop. 
This contrasts with more luminous radio sources (e.g. 
TN J1338) thought to be fuelled at higher rates through 
radiatively efficient standard accretion disks by cold gas 
(Best & Heckman 2012). These more luminous radio 
sources are hypothesized to have fuel brought in through 
mergers and interactions, which are in fact abundant 
in SPT2349—56. 'The debate thus remains open as to 
whether we are seeing a decaying radio source, or a radio 
source quickly building in luminosity. By better speci- 
fying the radio emission and its origin, we could learn 


about the build-up and state of the ICM that may al- 
ready be present at z — 4.3. 


5. CONCLUSIONS 


We have presented ATCA radio observations of 
SPT2349—56, a starbursting and gas-rich protoclus- 
ter, consisting of over 30 SMGs at z —4.3. We placed 
SPT2349—56 in context with uJy radio sources in the 
GOODS-N and ALESS fields, and with the other 22 
gravitationally-lensed SPT SMGs also observed with 
ATCA in our program. We also studied in detail the cen- 
tral galaxies identified by ALMA in SPT2349—56 near 
this strong radio detection. 

— We detected a single source at 2.2GHz in 
SPT2349—56, spatially coincident with the central three 
luminous members of the protocluster, denoted B, C, 
and G in Miller et al. (2018). While the ATCA radio 
centroid lies close to source C, which has the largest 
stellar mass in the protocluster, we cannot rule out that 
the radio emission is coming from B or G, or even a 
combination of the galaxies. 

- Under any of the possibilities above, the 214 uJy flux 
density at 2.2 GHz translates to more than 20 times the 
radio luminosity expected from the FIR-radio correla- 
tion defined by star-forming galaxies, and suggests that 
an AGN is driving the radio emission. 

- The radio source has a steep spectrum, with an in- 
dex of a= — 1.58 + 0.31, constrained by the ASKAP 
888MHz detection, and the non-detections at 5.5 and 
9 GHz, consistent with an AGN. 

- No other clear signs of AGN activity have yet been 
detected in this protocluster using any other diagnostics 
available to us (CO SLEDs; EW(OHies,um), [C11]/FIR 
ratios; optical spectra), highlighting the radio contin- 
uum as a powerful probe of obscured AGN in high-z 
protoclusters. 

- The three SMGs likely associated to the radio source 
have amongst the highest gas and dynamical mass of the 
protocluster members (Rotermund et al. 2021). More- 
over, high resolution ALMA imaging resolves this sys- 
tem into multiple interacting, star-forming clumps, with 
a surrounding arc of [C11] emission (Hill et al. 2020; 
Sulzanauer et al. in prep.). This is consistent with the 
idea that the availability of large amounts of gas and 
galaxy interactions, both of which are enhanced in gas- 
rich overdensities at high redshift, can trigger fast and 
obscured SMBH accretion. 

— No significant radio emission (nor any other robust 
AGN signature) is found from any other cluster member, 
constraining the radio-loud AGN content among SMGs 
in this protocluster to no more than 1096 (three of 30 
members), and likely just 3%. A radio stacking analy- 
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sis on the remaining ten brightest SPT2349—56 SMGs 
finds (11410) uJy, which is consistent with the average 
2.2 GHz emission from star formation via the FIR-radio 
correlation. We thus find no evidence that nuclear accre- 
tion powering radio emission exists below our detection 
threshold in other SMG members of SPT2349—56. How- 
ever, radio-loud AGN represent only 10% of all AGN, 
and X-ray observations and JWST infrared spectroscopy 
would be the next key steps to constrain AGN in this 
system and compare to AGN fractions found in other 
protoclusters. 

- The SPT2349—56 radio-loud AGN has a luminos- 
ity density of Lo 2 = 4.4 x 1025 W Hz-1, extrapolating to 
L14,rest = (2.4+0.3) x 102% W Hz^! with the measured 
a = —1.6, which is still over two orders of magnitude 
less luminous than the powerful radio galaxies normally 
studied at these redshifts. Many such HzRGs have rich 
protocluster environments, however it remains unclear if 
the opposite is true, that all massive z > 4 protoclusters 
have a central radio galaxy. 

- The fact that the radio AGN is detected in the hy- 
pothesized central seed of a growing BCG galaxy with 
significant stellar mass already in place makes this dis- 
covery an important new ingredient in understanding 
the formation and evolution of the cluster. 

- The radio luminosity was used to infer a radio jet 
power of Peay = (3.3 + 0.7) x 103° W, sufficiently large 
as to provide a dominant feedback on the cooling gas in 
the 10153 Mg halo. The radio luminosity also suggests a 
strong X-ray source with Lx = 10% W (integrated be- 
tween 2 and 10keV), easily detectable by Chandra or 
XMM-Newton. SPT2349—56 therefore has a high lumi- 
nosity AGN, even if in the form of a highly obscured 
quasar, and JWST will be a powerful tool to uncover 
its properties through high ionization infrared emission 
lines. 
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APPENDIX 


A. APPENDIX A 


In this appendix, we further assess the offsets between source centroids in the 888 MHz ASKAP image and the 
ATCA 2 GHz image that were discussed in section 3.1 (Fig 8). We measured peak fluxes in both images for all sources 
within a 13’ radius of SPT2349—56, measured their centroids, and calculated radial offsets for each. The offsets appear 
random in orientation, with the mean x and y offset being close to zero (0.3", —0.2"). In Fig. 8 we plot the radial 
offsets versus ASKAP flux. SPT2349—56 shows the largest offset, which could indicate that its origin may be physical. 
It has a 4.60 deviation from the median (excluding SPT2349—56) offset of 1.4". Even restricting the analysis to those 
sources with comparable flux densitites and SNRs (Ssss <2mJy) only increases the median offset to 1.6". 
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Figure 8. Assessing the offsets between ASKAP and ATCA sources. Top: The 888MHz ASKAP image surrounding 
SPT2349—56 (green LABOCA contours, showing ALMA sources as green dots) with ATCA 2 GHz contours overlaid (26' x 19' 
field shown). Bottom: The radial offsets between the centroids for sources common in both images. SPT2349—56 shows 
the largest offset, which might therefore require a physical interpretation. It has a 4.60 deviation from the median (excluding 


SPT2349—56) offset of 1.4" offset. 
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B. APPENDIX B 


In this appendix and Figs. 9 — 11, we show the CO J — 7, 11, and 16 lines for the central B, C, and G sources, whose 
line strengths are plotted in the SLED diagram (Fig. 5). We also show the H20 lines that are used to compare with 
the FIR luminosity estimates from Hill et al. (2020). 

In order to measure line strengths, the bright and well-detected [C 11] lines provided in Hill et al. (2020) were used 
as a template. These [C r1] lines were fit by single and double Gaussian profiles, and we selected the integration range 
by scaling the [C ri] profile to the rest frequency of the line of interest and then summing channels between —2o and 
20 (where o is the standard deviation of the best-fitting linewidth), or for cases where two Gaussians were a better 
fit, from —2o1, to +20R, where o1, and op are from the left and right Gaussian fits, respectively. 

The CO(7-6) line is blended with the [C1](2-1) line, and the CO(16-15) line is blended with the OH doublet, so 
these had to be fit and subtracted before integrating over the CO lines. For the former case, where both CO(7-6) and 
[C1](2-1) are both well-detected in B and C, we simultaneously fit single Gaussian profiles at the locations of the two 
lines, then subtract the best-fit [C1](2-1) model from the spectrum, and sum over the relevant channels as described 
above (then vice-versa to obtain [C 1](2-1) line strengths). For source G, we do not see any strong line features around 
the expected [C1](2-1) frequency, so we simply sum over the CO(7-6) and [C 1](2-1) channels in the raw spectrum. 
The CO(16-15) line is not well-detected for any sources but the OH doublet is, so we fit a Gaussian to these OH 
lines and subtract the models before summing over the CO(16-15) channels. In the fit we force the amplitude of each 
doublet component to be equal, and we fix the width of each doublet component to be equal to the width of the [C 11] 
line (described in Sec. 3.3; see Table 4). Since the profile for G is two Gaussians, we include an additional OH doublet 
component of equal amplitude and fixed frequency separation/width to match the [C r1] profile. This leaves two free 
parameters in all fits: the frequency of the first doublet, and the amplitude of the all the components. In Fig. 9 we 
can see that for source B the two OH doublet components are blended with each other due to the large FWHM of the 
system, and for G the four components blend into three peaks. 

Lastly, Band 4 and Band 6 continuum flux densities were estimated by averaging over all line-free channels in the 
original (non continuum-subtracted) data cubes (again using the [C 11] line as a template). We combined channels from 
the lower and upper sideband of these observations, meaning they are at observed frequencies of 147 and 231 GHz, 
respectively. Band 7 continuum flux densities (around the CO(16-15) and OH lines) are already provided in Hill et al. 
(2020). 


C. APPENDIX C 


Here we describe the ATCA observations of the full sample of 23 SPT-SMGs observed in the survey program, shown 
in Fig. 3. These SPT-SMGs were drawn from the complete sample of 81 sources (Reuter et al. 2020), selecting those 
that had the best redshift constraints at the time of observations. All but three of the 23 SPT-SMGs are detected 
at > 4c significance at 2.2 GHz. The 2.2 GHz flux densities are measured at peak pixels (Table 5), as in all cases the 
sources are unresolved in the 8" x 5" beam. The restored beam sizes and position angles are also listed in Table 5. 
ALMA 850 um overlays are shown in Fig. 12 (data from Spilker et al. 2016; Reuter et al. 2020). LABOCA 850 um 
fluxes and ALMA-derived redshifts from Reuter et al. (2020) are also listed in Table 5 for completeness. 

Most sources are not detected or only marginally detected at 5.5GHz (nine detections at > 30) and 9.0 GHz (four 
detections at > 30). For those sources detected at these higher frequencies with ATCA, we measure flux densities from 
peak pixels when the source is unresolved, or as aperture measurements when the source is resolved. We show the 
nine sources detected at 5.5 GHz in Fig. 13, the four sources detected at 9.0 GHz in Fig. 14. We have also searched for 
detections in the ASKAP 0.888 MHz RACS survey described in section 2.2, listing their flux densities in Table 5. We 
find 15 of the 23 sources are significantly detected by ASKAP. 

We derive radio spectral indices directly for all sources with at least two radio detections, and list these together with 
flux densities in Table 5. The data was fit according to a linear function using a Markov Chain Monte Carlo algorithm 
(MCMC) implemented by the emcee package (Foreman-Mackey et al. 2013). This MCMC package samples the posterior 
probability function, and is used to determine the error contours shown in Fig. 2, as well as the uncertainties on a in 
Table 5. We summarize their radio spectral indices in in Fig. 15 and Table 5. 

'The lensed SMGs are shown in Fig. 3, where we estimate their rest 1.4 GHz luminosities directly using the measured 
a, or with a = — 0.8 if only detected at a single radio frequency. In general these SPT-SMGs follow the same FIR-radio 
correlation as the other field samples shown. However, three sources are highly significant outliers from the FIR-radio 
correlation: SPT0125—50 at z=3.96; SPT0202—61 at z=5.02; and SPT0550—53 at z —3.13. Given how rare such 
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strong outliers are in the field SMG samples (only one of 76 SMGs shows anywhere near this level of radio excess 
in the ALESS SMG sample - Thomson et al. 2014), we propose that the lensing galaxy rather than SMG may be 
the more likely radio-AGN in these three cases. These radio excess sources exhibit steeper radio indices than typical 
star-forming galaxies, comparable to or exceeding SPT2349—56. Without knowing if the lens or source redshift is 
correct, we cannot reasonably apply the radio K-correction to estimate the rest 1.4 GHz luminosity, and therefore we 
do not include these three in figure 3. 

In particular, SPT0550—53 shows an extended radio morphology /jet, well resolved in all three ATCA frequencies, 
which is more naturally explained by a lower redshift radio-loud galaxy. Further, the optical spectrum of the lens 
SPT0550—53 shows AGN emission lines. Neither SPT0125—50 nor SPT0202—62 show AGN signatures in their optical 
spectra. SPT0125—50 is curious as the ATCA 2GHz flux density is very close to that expected from the FIR- 
radio relation, however ASKAP reveals a 3mJy source well centered on SPT0125—50, implying an incredibly steep 
a = —2.24. Thus any K-correction to lower rest frame frequencies than that probed by 2.2 GHz observations quickly 
places SPT0125—50 significantly above the FIR-radio correlation. 

These results also beg the question of whether the radio emission in other lensed SPT-SMGs might be contaminated 
from the often massive lens galaxy (Rotermund 2020). While in some examples, especially in SPT0538—50, the radio 
emission is directly identified as coming from the ALMA-detected lensed SMG components, in others the Einstein 
radius of the lensed source (Spilker et al. 2016) is too small to be detected offset from the lens galaxy itself, even 
at 9GHz. The distribution in a constrained by the fits in figure 15 show a mean of —0.93 + 0.14, offset steeper, 
but still consistent, with the a = —0.8 found in samples of unlensed SMGs reported in section 3.2 (e.g., Thomson 
et al. 2014). Several of the higher redshift sources in figure 3 do in fact show a marginal excess over that expected 
from the FIR-radio correlation. This excess sometimes appears only due to the comparison shown at rest 1.4 GHz, 
accentuating the K-correction from their steeper than average a we measure. Given the large uncertainties from the 
often two-point a estimates, this may be inconsequential. Generally, optically identified AGN are relatively rare in 
the distant red galaxies that are often found to lens these SMGs (Rotermund 2020), and the gravitationally boosted 
radio signal associated with the high-SFR SMG is a more probable source of the strong radio emission we see in these 
19 SPT-SMGs. Their radio emission is not obviously contaminated by their foreground lens. 
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Figure 9. Cutouts and spectra of CO(7-6), CO(11-10), CO(16-15), and H20 line emission for galaxy B. The cutouts in each 
panel show continuum emission (obtained by averaging over all line-free channels) and line emission (obtained by averaging over 
all channels where the line is expected — see Section 2 for details), with contours starting at 20 and increasing in steps of 3c. 
Apertures are shown as red circles, and used to obtain the spectra shown in the right panels. In each spectrum plot, we show 
the [C11] profiles from Hill et al. (2020), scaled to the expected frequency of the given line, and arbitrarily normalized. The 


shaded regions show the integration ranges (set to be 4 


-2e about the [C1] line — see Section 2) used to obtain line strengths. 


The CO(7-6) line is blended with the [CI] (2-1) line, and the expected central frequency (or for G, two central frequencies as the 
[C 11] profile has two components) of the [CI](2-1) is marked in red. The [CT](2-1) line is fit by a Gaussian profile and subtracted, 
and the original spectra are shown by the dashed lines. Similarly, the CO(16-15) line is blended with the OH doublet, and we 
mark the mean frequency of each OH line in red (corresponsing to two frequencies for B and C, and four frequencies for G). 
The OH doublet is fit by a scaled [C 11] profile (see Section 2) and subtracted, and the original spectra are shown by the dashed 


lines. 
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Figure 10. Same as Fig. 9 but for galaxy C. 
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Table 5. ATCA observations at 2.2, 5.5, and 9.0GHz, as well as ASKAP at 888 MHz, of 23 SPT SMGs. We also include 
unresolved LABOCA 850 um flux densities (S850), ALMA redshifts, and best-fit spectral indices. 


Name tint S2.2 Beam" PA? Si. ium Si, S850 z a Comment! 
(hrs) | (uJy) | (77) (deg) (uJy) —(uJy) (uJy) (my) 

SPT0027—50 0.58 334431 8x4 106 «105 «150 1064201 48 3.444 —1.28+0.25 nn 

SPT 0103—45 0.54 229429 8x4 16 «111 «159 « 589 125 3.092 - nn 


SPT0109—47 0.56 1106-231 9x5 11 824436 461452 1417+194 109 3.614 —0.42+0.04 yR yR 
SPT 0125—47 0.59 586435 9x5 112 382236 193+49 1835+201 144 2.515 —0.82+0.10 ymR 


SPT 0125—50 0.56 365429 5x7 110 <102 <156 2792+194 109 3.959 —2.25+0.12 nn 
SPT 0202—61 0.68 710435 9x5 6 22540 <153 19434179 109 5.018 —1.16+0.09 yn 
SPT 0245—63 0.71 9433 9x4 TT <99 <135 < 598 61 5.626 = nn 
SPT 0345—47 0.68 275429 8x6 174 <99 <132 10143194 89 4.296 —1.03+0.38 nn 
SPT0346—52 0.72 1621238 5x9 68 «105 «132 « 603 131 . 5.656 = nn 
SPT 0418—47 0.70 173422 8x5 163 <93 <135 430198 108 4.224 —0.99+0.67 nn 
SPT 0512—59 0.56 465434 7x5 153 <177 <231 1531-197 75 2.233 —1.32+0.17 nn 
SPT 0529-54 0.54 260450 8x5 153 <120 <180 < 586 118 3.369 B nn 
SPT 0532—50 0.57 489439 8x5 154 144449 «177 109340231 118 3.399 —1.06+0.16 mn 


SPT 0538—50 0.56 581536 8x5 157 341-59 168458 1490+237 125 2.786 —0.81+0.13 yR mR 
SPT 0550—53 0.55 1288448 9x6 169 446439 270456 4060+187 53 3.128 —1.22+0.05 yR yR 


SPT 0551—50 0.56 286425 8x6 160 159445 <171 5201-185 74 3.164 —0.70+0.24 mn 
SPT 2031—51 0.48 269431 9x5 51 <123 <186 7214203 65 2.452 —1.08+0.39 nn 
SPT 2134-50 0.53 334447 8x4 33 174443  <162 8043-196 101 2.780 —0.90-£0.44 yn 
SPT2319—55 0.54 75444 8x4 4T «126 «174 « 600 38 5.293 = nn 
SPT 2332—53 0.56 244423 9x5 36 146441 <162 < 581 57 2.756 —0.82+0.56 mn 
SPT 2349-56 0.56 215227 8x4 30 <120 <162 867189 1067 4.303 —1.580.31 nn 
SPT 2353—50 0.56 24453 9x5 30 <138 <159 < 594 A1 5.576 -= nn 
SPT 2357—51 0.55 131419 9x5 31 <108 <156 < 589 53 3.070 = nn 


The 2.2 GHz beam is quoted as x and y FHWM. The 5.5 GHz beam is typically 3.6” x 2.2”. The 9.0 GHz beam is typically 
2.2" x 1.2”. 

^ The PA of the 2.2 GHz beam is the angle East of North. 

“ The radio spectral index a, defined as S, «v^. 

d In SPT2349—56 we have assumed source C with Sgso um — 4.7 mJy is the host of the AGN, although it could be B or G as 
described in the text; here we still provide the unresolved LABOCA flux density. 

* At 5.5 and 9.0 GHz, the 3c upper limit is listed unless there is a detection at > 3o at the ALMA position. 

îi The 888 MHz measurements are from the ASKAP RACS survey, described in section 2.2. Sources with < 30 positive signal 
are listed at these limits. 

tit Comments list whether 5.5 GHz and 9.0 GHz data show detections > 4c (y), marginal detections (m) where < 4o flux 
density is measured at the ALMA position, or no detection (n). We indicate the four sources with resolved radio morphologies 
(R), in SPT0538—50 clearly following the ALMA emission, although in SPT0109—47 and especially SPT0550—53, the resolved 
emission appears to be an extended lobe or jet. 
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Figure 12. ATCA 2.2 GHz images (30” x 30”) of the 23 SPT-SMGs observed, with ALMA 850 um contours overlaid. 
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Figure 13. ATCA 5.5 GHz images (30” x 30”) of the nine detected SPT-SMGs, with ALMA 850 um contours overlaid. 
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Figure 14. ATCA 9.0 GHz images (30" x 30") of the four detected SPT-SMGs, with ALMA 850 um contours overlaid. 
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Figure 15. Radio spectral indices with errors constrained by MCMC modelling for gravitatinally lensed SPT-SMGs having 


at least two detections between the ATCA and ASKAP followup. Grey shaded regions show the 1 and 20 uncertainties on a 
derived from the ATCA data. The brighter sources with steeper spectra are generally detected by the ASKAP RACS survey. 


The fit for SPT2349—56 is shown in Figure 2. 


